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We study the state decay of two qubits interacted with a common harmonic oscillator reservoir. 
There are both decoherence error and the error caused by the amplitude change of the superradiant 
state. We show that frequent 7r-phase pulses can eliminate both typpes of errors therefore protect 
a two-qubit odd-parity state more effectively than the frequent measurement method. This shows 
that the the methods using dynamical decoupling and the quantum Zeno effects actually can give 
rather different results when the operation frequency is finite. 
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I. INTRODUCTION 

The interaction between a quantum system and its en- 
vironment inevitably leads to the decoherence [l|, Q of 
a quantum state. Such quantum decoherence can of- 
ten cause severe distortion to a quantum state render- 
ing many quantum systems in the real world usclcss[3- 
Il0| . In order to protect a quantum state, many methods 
against decoherence have been studied. Among the ex- 
isting proposals, most of them are for single-qubit state 
protection [ll| - fl6| . Recently 17], a scheme for the protec- 
tion of quantum entanglement of two qubits at OK tem- 
perature was proposed using quantum Zeno effect (QZE), 
i.e., via frequent measurement of the environment photon 
number for the Jaynes-Cummings (J-C) model. However, 
as shown below, besides the decoherence error, the am- 
plitude of the superradiant state can also changed. It 
decreases with evolution. Intuitively speaking, the su- 
perradiant state changes into |00) gradually in the evolu- 
tion therefore the initial odd-parity state cab be severely 
distorted after a long time evolution. The frequent mea- 
surement method cannot eliminate such distortion effi- 
ciently because it actually removes the term 1 00) at every 
step. After a long evolution time with a fixed measure- 
ment frequency, the amplitude of the superradiant state 
decreases a lot therefore severely distort the initial un- 
known state. On the other hand, given the existing tech- 
nologies, it seems that the measurement of photon num- 
ber of the environment of all modes remains a challeng- 
ing task. It is therefore an interesting problem to study 
how one could protect a two-qubit state with techniques 
which have been demonstrated already, for example, dy- 
namical decoupling scheme [llTfl6l fl8f which have been 
demonstrated experimentally rece ntly[llll|. It is well 



known that in the limit of infinitely frequent operations, 
QZE and dynamical decoupling are unified and can have 
the same results [22rl24j |. The two methods are not com- 
pared in the more realistic condition when the operation 
frequencies are finite. Here we show that the dynami- 
cal decoupling scheme achieved through a frequent ap- 
plication of 7r-phase pulses can protect a two-qubit state 
more effectively. The scheme not only protects the state 
from decoherence error, but also prevents the amplitude 
changing of superradiant state. 

This paper is arranged as following: we first review the 
existing results of the J-C model[l7|,[29| for two qubits, in 
particular, the time evolution of the odd parity state |17j 
under zero temperature. We point out why the ampli- 
tude of superradiant state changes in a frequent mea- 
surement scheme. We then show how to protect the state 
by 7T— phase pulses and why n— phase pulses scheme can 
prevent the amplitude change. The consequences of the 
finite frequency and duration of each pulses are also pre- 
sented. 



II. AMPLITUDE CHANGE OF 
SUPPERRADIANT STATE IN J-C MODEL 

Consider the following Hamiltonian for a two-qubit 
system and its environment as used in}17j|: 



H — H s + H e + Hi 

where 

H s = LJ (crf(Ti + o\a^) 
is Hamiltonian of the two-qubits (system), 
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is the Hamiltonian of the environment, and 
Hi = (aicr^ + cr 2 cr^) ^ gkbk + h. 
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is Hamiltonian of the interaction between the system and 
the environment. Notations 6^ and are the annihila- 
tion and creation operators of the environment with fre- 
quency Uk ; loq is the atomic transition frequency between 
the ground state |0) and the excited state |1); a + = |1)(0| 
and a~ = |0)(1|. The solution of such a model for the 
case of zero environmental t emp erature is well known 
and it can be found in Ref . [29|, |3fJ ■ In particular, for 
an odd-parity two-qubit initial state, there exists a dark 
state 

|M> = -|1)i|0>2--|0} 1 |1> 2 (5) 
a a 

and a superradiant state [25l - [28| . 

|i/) = -|l)i|Q) 2 + -|0)i|l) 2 . (6) 
a a 

The dark state \fi) does not change with time under the 
JC model, while the superradiant state changes with time 
according to 

k)®|0)e=T?(i)|^)®|0)e + |0) 1 |0) 2 ®^(c(t) fe |l fe ) e ) (7) 

k 

and 

7j(t) = e-r [cosh(— ) + - sinh(— )] (8) 

where f2 = \A 2 - 4i? 2 , a = y/af + a\ and R = aW. 

Given any odd-parity initial state \ipo) = + /3 2 |f), 
the time evolution is 

Pi(t)=Pi (9) 

fo(t) = fav(t) (10) 

and rj{t) is given in Eq.©. 

\ip(t)) = ( ( 8i|M)+ / 8 2 r?(*)k»<8|0>e+^l )il°)2®c(*)fc|lfc}e 

k 

(11) 

The decoherence error come from the term ^ fc |0)i|0) 2 ® 
c(t)fc|lfe) e - As shown in Ref. 17], by frequently measur- 
ing the environment, one can remove the the term |00) 
and protect the initial state from decoherence error, as 
long as one does not find a photon coming from the reser- 
voir. Intuitively speaking, such a frequent measurement 
works like a state filter which removes 1 00) during the 
stage when its probability is small. However, even though 
one can always remove the term with 1 00) successfully by 
measurement, one can not protect the initial state for a 
long time with the scheme because r/(t) decreases signifi- 
cantly with time. Suppose the environment is measured 



after every time interval At, and we continue to find no 
photon. At time i, the state is changed into 

\m) = (/3i|//>+/3 2 r(*)W)®|0)e (12) 

and 

r(t) = [v(At)] t/At . (13) 
Since each measurement removes the photon in environ- 
ment, |f) <S> |0) e restart the evolutionm from the initial 
state again in each time interval. To protect two-qubit 
state of the system more effectively, we can use the dy- 
namical decoupling scheme through applying tt pulses 
frequently instead of a filtration scheme with frequent 
measurement. The main idea is, whenever the state de- 
cays a little bit, i.e., r](t) decreases a little bit and the 
term with 1 00) appears with a small amplitude, a n pulse 
is applied and as a result, rj (t) will rise and the amplitude 
of term with |00) decreases. That is to say, a tt— pulse 
does not simply remove the term with |00), it changes 
the term with 1 00) back to the superradiant state. There- 
fore, it differs from the state filtration of the measurement 
based scheme - it is really a more effective scheme of state 
recovery. 



III. ELIMINATION OF DECOHERENCE WITH 
FREQUENT tt-PHASE PULSES 

We show here that 7r-phase pulses can eliminate the de- 
coherence on the one hand and prevent amplitude chang- 
ing of the superradiant state on the other hand. A tt- 
phase impulse take a phase-shift operation as: 

|0> ^ -|0> |1> -> |1> (14) 

Apply a 7T— pulse to each qubit, the two-qubit unitary 
operation is then 

|l)i|0) 2 ->|l>i|0) a 

|0)i|l) 2 -»• |0)i|l) a (15) 
|0)i|0) a -> -|0>i|0) a 



A. Some iteration formulas 

Our method involves applying simultaneously 7r-pulses 
to each qubit. To obtain the results of the method, we 
need some iteration formulas first. Suppose the interval 
between two consequent impulses is r and n = [^], we 
can calculate the coefficients, r\ and r 2 , by 
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pz ■' />(m+l)r 

(*) = - / dnf(t - «)[a?ri(») + aia 2 r 2 ( K )] - ^ (-1)"-™ / d«/(t - «)[a?ri(fe) + aia 2 r 2 (/s)] 

Jut m _ n 7mr 



m=0 
n-1 



r 2 (i) = - / <f/s/(t - «)[alra(«) + aia 2 ri(«)] - ^ (-1)™-™ / cte/(f - ^[^(k) + aia 2 ri(/c)] 



(16) 



where we know c(t)k, the coefficient of |lfc) e in Eqn.|llj. 
changes the signal because of the two phase impulses cor- 
responding to t = mr(m is integer.). 

Under the interaction with a Lorentzian spectral den- 
sity, we get the similar result as the free evolution, namely 
that there exist one dark state \fi) and one superradiant 
state \v). We next study the condition that fi\(t) and 
/3 2 (t) should satisfy between one time interval. 

When t is between mr and (m + l)r. 

0i(t)=A (17) 
P 2 {t) + \p 2 {t) + R 2 (3 2 {t)=0 (18) 
We know that the general solution of Eqn. [TBI is: 



'[An COSh(- 

y Q(t — mr) 
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)] 



(19) 



, sinh(- 

where t € [mr, (m + l)r]. A m and B m are the constant 
coefficients of the solution f3 2 (t) at each time interval. 
Our task now is to determine the relation of all A m and 

- 

At t = (m + l)r using the boundary condition, (3 2 ((m+ 

1)t~) = A((m + l)r+) and ft((m + l)r-) = -/3 2 ((m + 
1)t + ), we can get the relation: 



Ai 



[A 



.cosh( — ) 
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B TO sinh(— )] 



-l = - e 



vr fir fir 

2 [A m sinh( — ) + B m cosh( — )] 



(20) 



— e ^ [A ro cosh(— ) 



• smh (^r)] 



Here j3 2 {mT~) means the value of /3 2 (mr) before the 
pulse. And /3 2 (mT + ) is the value of /3 2 (itit) after the 
pulse. 

We know that if the initial state is superradiant state, 
(3 2 (t) — 1, the initial value should be A — 1 and -Bn = A 
By using the relationship between A m 
B m +i, we can get an analytical function £(t) which can 
express the fidelity 
state. 

e -^^[^ m cosh(^l)) 



and 



\J (v\p v {t)\v) of superradiant 



-B m sinh( 



— mr) 



)] 



(21) 



for t € [mr, (to 
is given by 



l)r]. The fidelity F(t) = ^ (yj\p(t)\yj) 



F(t) 



l/3i| 2 + l/3 2 | 2 ^) 



(22) 
(23) 




0.985 

F(t) 

0.98 



0.965, 



t 0.5 



FIG. 1: (Color online) Comparison of the fidelity evolution 
during t = [0, 1] of two methods for the superradiant state. 
We set r = 0.1, A = 2 and fi = 1. Line 2 is for dynamical 
decoupling of this work and Line 1 is for QZE of Ref[l7(- The 
fidelity only oscillates in a small range around a horizontal 
line in the dynamical decoupling method. 



With this iteration formula for £ (t) above, we see tah the 
fidelity oscillates about a horizontal line within a small 
range after a pulse is applied, first increasing and then 
descending, as shown in Fig pQ. In contrast, the fidelity 
based on the method of Quantum Zeno Effect descends 
almost monotonously with time. 



Effect of finite duration of Double-7r-Phase 
Operation 



In practice, one cannot set the duration pulse time to 
be infinitely small. Here we consider the more realistic 
case that duration time of double- 7r-phase is finite and we 
study show the effect on the fidelity. We only consider 
the sequential pulses method and we suppose the pulse 
duration is in each time interval r. We know that the 
decoherence coefficient of Superradiant state is Eqn. [TT)1 
when t G [tot, (to+1 — t^)t]. So the effective Hamiltonian 
for the sequential double- 7r-phase operator is H p hase — 
^(a+af + <7 2 + a 2 -)E ro > 1 [e(mr - - 9(mr)] and 
@(x) is the step function. Under the Hamiltonian H a u = 
H s + H e + Hi + Hphase, we can write down the state 
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during the time of [(m 



jr)r, tut]: 

iJ ^ (t - (ro -* )T) (r(t)i|l>i|0)i 
r(t)2|0)i|l) 2 )®|Q) e 
X)|0>i|0> 2 ®c(t) fc |l fc ) e 



more effective in protecting the states. Our result clearly 
shows that the quantum Zeno effect and the dynamical 



(24) 



In the same way as shown above, we also can get a inte- 
gral equation describing the coefficient r\ = — J * drf'it— 
T)\a\ri(T) + aiot2r<2(T)\ and the new correlation function 
of the duration time should be f'(t) — W 2 e ( -^ x+tE ^ L ^ t . 
Using the boundary condition of t = (m — -h)T and 
t = tot, we eliminate the coefficient in the operator's du- 
ration time and get coupled relationships between A m , 
B m and A m+X , B m+ i(Eqn[T7 
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FIG. 2: (Color online) The initial state is the superradiant 
state. Line 2 is the fidelity under sequential double-7r-phase 
)t)1 method (Dynamical Decoupling). Line 3 is the fidelity under 
QZE. Line 1 is the fidelity of free-evolution under the inter- 
action between system and environment. From the result, we 
can easily find the sequential-pulse method is much better 
than the QZE. 
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where 7 = 

The time evolution of the system is given by \ip(t)) 
(/3i|/*>+/82«t)k»®|0)e+E*|0)i|0) 2 ®c(t) fc |l fc ) e aiid€(t) 
is defined by Eqn. 19 with the new relation of Am and B m 
above. Also, the fidelity of the state under decoherence 
with the original state is 



F{t) = \^ + \(3 2 \^{t) 



(26) 
(27) 
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C. Numerical Calculation of Dynamical 
Decoupling and Free-evolution 

In case 1, we set r = 0.1, A = 2 and fi = 1 and 
we compare the fidelity of the free-evolution, sequential 
double- 7r-phase operator and random one. The result is 
shown in Fig[5]. In case 2, we set t = 0.2, A = 2, Q. = 1, 
Ni = 10 and N2 = 20 to compare how the duration time 
of pulses affect the fidelity. The result is shown in Fig [3], 



IV. CONCLUDING REMARK 

We present a strategy to protect the odd-parity states 
of two qubits under temperature environment by fre- 
quently applying the 7r— pulses. Comparison between this 
method and the method based on frequent measurement 
is done, it seems that the frequent- tt— pulses method is 



FIG. 3: (Color online) The initial state is the superradiant 
state. Line 2 is the method with the instantaneous pulse. 
Line 4 is for N = 20 and line 3 is for N = 10. Line 1 is for 
the free-evolution. 



decoupling can have rather different results when the op- 
eration frequency is finite, though the two methods give 
essentially the same results in the limit of infinite opera- 
tion frequent as shown in Ref. [22h24| . 
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